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ABSTRACT 

We have previously presented the results of a 325 MHz radio survey of M31, conducted with the 
A-configuration of the Very Large Array. In this survey, a total of 405 radio sources between <6" and 
170" in extent were mapped with a resolution of 6" and a sensitivity of ~0.6 mJy beam~^. Here, we 
compare the resultant source list and image with other radio, IR, optical, X-ray observations and catalogs 
of the region. Through this, we were able to identify five supernova remnant (SNR) candidates and three 
pulsar wind nebula (PWN) candidates in M31, as well as three Milky Way radio stars, a possible Milky 
Way Planetary Nebula, and a bevy of interesting extragalactic objects: a BL Lac, a Giant Radio Galaxy, 
a galaxy merger, and several high-z radio galaxy candidates. In addition, a large number of compact 
(9 < 6") extremely steep spectrum sources (a < —1.6; Si, oc v°') were detected whose nature is unknown. 



Subject headings: catalogs — galaxies: individual (M31) 
continuum: galaxies 



radio continuum: general — radio 



1. INTRODUCTION 

The radio population of a star-forming galaxy - dom- 
inated by supernova remnants (SNRs), Hii regions, and 
radio pulsars - contains much information on the evolu- 
tion and properties of such galaxies. To better understand 
such sources in M31, we have surveyed M31 at 325 MHz 
with the Very Large Array^ (VLA) in order to identify 
radio sources intrinsic to this galaxy. In the introduction 
of ? ? (hereafter referred to as Paper I), we discussed 
some of the advantages of surveying M31 at low frequen- 
cies, mainly that the large field of view of the VLA at this 
frequency allows one to image the optical disk of M31 in 
one pointing. However, observing at low frequencies also 
has advantages in classifying the detected sources. Low 
frequency observations (z^ < 1 GHz) have increased sensi- 
tivity to the steep-spectrum (a < —0.5) population of a 
galaxy, mostly composed of SNRs and radio pulsars, than 
higher frequency observations. These two different class of 
objects provide separate pieces of information: the SNR 
population traces recent massive star formation and indi- 
vidual SNRs both highlight the local interstellar medium 
(ISM) and are crucial to understanding the cosmic ray 
population of the galaxy, while radio pulsars are very use- 

^ The VLA is operated by The National Radio Astronomy Observatory, which is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 



ful in determining the line-of-sight electron density and 
magnetic field and contain information on the total core- 
collapse supernova rate of a galaxy. In a low-frequency 
survey, these sources are noticeably brighter than the flat- 
spectrum population of a galaxy, composed mainly of HII 
regions and Pulsar Wind Nebulae (P WNe) . By comparing 
our images to higher frequency images of M31, we can sep- 
arate these different components of the radio population 
of M31 from each other. 

In addition, low-frequency observations coupled with 
higher-frequency ones allow the detection of spectral cur- 
vature, a powerful probe of a source and its surround- 
ings. Two processes which can cause low-frequency spec- 
tral turnover are free-free absorption, indicative of ion- 
ized material along the line of sight, and synchrotron self- 
absorption, which is believed to occur in active galactic 
nuclei (AGN; ? ?). Simultaneous 74 and 333 MHz obser- 
vations have detected free-free absorption caused by ion- 
ized gas inside SNR Gas A (?), and 74 and 327 MHz obser- 
vations of SNR W49B discovered a foreground HII region 
obscured by a molecular cloud (?). Low-frequency surveys 
of the center of star-burst galaxy M82 allowed the authors 
to calculate the properties of its ionized ISM through the 
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detection of free-free absorption (e.g. ? ?; ? ?). In this 
case, unhke for the the Milky Way cases discussed above, 
spectral turnover was detected above 325 MHz. Low fre- 
quency spectral turnover has also been observed in com- 
pact, distant, powerful radio galaxies known as Gigahertz 
Peaked sources (GPS) believed to be examples of AGN 
powered radio galaxies whose jets have yet to penetrate 
through the ISM of the galaxy (?). In addition, a radio 
spectrum with low-frequency values above the power-law 
extrapolation of higher frequency nicasurcnients is the sig- 
nature of nonlinear shock acceleration, an important pro- 
cess in explaining the observed cosmic ray spectrum, and 
has been observed in young SNRs such as Tycho's SNR 
and Kepler's SNR (?). 

In Paper I, wo described a 325 MHz ("GLG") survey of 
M31 conducted with the A-configuration of the VLA. In 
this survey, we imaged 405 radio sources between 6" (the 
beam size) and 170" in extent, and calculated for each 
source its angular distance from the center of M31 {Om'si ), 
the projected radius of the source assuming it is in the 
optical disk of M31 (Rmsi), its major axis {9m), minor 
axis {9m), position angle {9pa), peak intensity {A), inte- 
grated flux {S), and - through comparisons with the 36W 
(?), 37W (?), Braun (?), and NVSS (?) radio surveys - 
its spectral index a and spectral curvature parameter (p, 
using the convention: 

5^ oc z."e-'^('^"'"') (1) 

To determine if a "GLG" source had a counterpart in one 
of these other catalogs, we used a limiting distance of: 



Max 



oGLG I oGLG 



(2) 



where 6'^'^ is the major axis of the GLG source, ^J^^*^ is 
the minor axis of the GLG source (these are not the decon- 
volved values but the size of the source in our images) , 9m 
is the deconvolved major axis of the catalog source, and 
9m is the deconvolved minor axis of the catalog source 
(?). Only a "good match" - one where the GLG source 
had only one match in the other catalog and its counter- 
part only had no other matches in the GLG source list - 
was used to calculate a and f. 

In addition, we classified each source into one of four 
morphological categories: 

• Unresolved ( "U" ): The 281 GLG sources with Om < 
2 X 9m and 9m < 15". 

• Elongated ("El"): The 16 sources with 9m > '2 x9m 
and 9m < 15". 

• Complex ("C"): The 51 sources with 9m > 15" that 
have structure on scales smaller than ~10". 



one of the aforementioned higher frequency radio 
surveys. Extended sources without counterparts in 
one of these surveys are most likely spurious, but 
were included in the final source list for complete- 
ness. 

This source list is providc;d in Tabic 3 of Paper I. 

As described in Paper I, a statistical analysis of the an- 
gular distribution (^m3i) of sources, the radial distribution 
of sources {Rm^h), flux density (S'325), and spectral index 
{a) distribution of sources revealed that the majority of 
detected sources are background radio galaxies. The pri- 
mary purpose of this paper is to identify sources intrinsic 
to M31, determine their nature, and discuss their proper- 
ties. In addition, we have also identified interesting back- 
ground and foreground sources. This paper is structured 
as follows: Section 2 discusses the information used to 
classify sources. Section 3 presents the final classification 
of sources, and Section 4 summarizes our results. 

In this paper, 5*,^ is the flux density of a source at a fre- 
quency v MHz, and the distance to M31 is assumed to be 
780±13 kpc (?). 

2. CLASSIFICATION OF OBSERVED SOURCES 

In Paper I, we provided statistical information which 
helped determine the general properties of the GLG 
sources. However, we also want to determine the natiirc of 
each observed source in order to sort out background radio 
galaxies from sources in M31 and the Milky Way, as well as 
to differentiate between the different kinds of sources that 
comprise these categories. In order to accomplish this, we 
used three different pieces of information: a source's ra- 
dio spectral properties, its morphological characteristics, 
and its multi-wavelength (from the far-IR to the X-ray) 
properties. While it was not possible to determine the 
nature of all 405 GLG sources, we were able to classify 
125 of them based on the information and methodology 
described below. 

2.1. Radio Spectrum Classifications 



• Extended ("Ex"): The 57 sources with 6'm > 
15" that do not have structure on scales smaller 
than ^10". These sources were detected by running 
the source-finding algorithm SFIND, part of the 
MIRIAD software package, on images smoothed to 
20"x20". A subset of these category are the Good 
Extended Sources ("ExG"), the two sources which 
satisfy the above criteria and have a counterpart in 

2 Only "U", "El", "C", and "ExG" 
of Paper I for more details. 



We defined three spectral categories in order to help 
separate different classes of sources steep spectrum {a < 
—1.3 and \ ^ \< 3a ^, where a^^ is the error in tp), flat spec- 
trum (—0.5 < a < 0.5 and | <p \< 3(7,^), and sources with 
spectral curvature ( | </? |> 3(t^).^ In the cases of sources 
with spectral curvature, the source must have been de- 
tected at two other radio frequencies - most commonly 
610 and 1400 MHz. 105 GLG sources fell into at least 
one of these categories. Table ?? lists the steep spectrum 
sources. Table ?? lists the flat spectrum sources, and Ta- 
ble ?? lists the sources with spectral curvature. Figure 1 
shows the positions of these sources with respect to the 
optical disk of M31. The distribution of steep spectrum 
and flat spectrum sources is uniform on the sky, while the 
distribution of turnover sources is condensed around the 
optical disk of M31. However, this most likely is a selec- 
tion effect - as mentioned in Paper I, (p was only calculated 
for those GLG sources with counterparts in at least three 
other radio surveys. As shown in Figure 1 of Paper I, these 
surveys only overlap around the optical disk of M31. 

sources are included in these categories, since the 55 "Ex" are unlikely to be real sources. See Section 2.3 
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This classification scheme was defined in such a way to 
separate different physical classes of sources. 

Steep Spectrum Sources With the advent of deep 
low-frequency surveys, Ultra-Steep Spectrum sources 
(USS; a < —1.3) have been increasingly studied and were 
discovered to be efficient tracers of High-z Radio Galaxies 
(? ?; ? ?) and relic radio galaxies (? ?; ? ?). In addi- 
tion, pulsars - cither in M31 or in the Milky Way - would 
fall in this category, since pulsars have an average spectral 
index of a ^ —1.8 (?). Thirty-six non-Extended sources 
fall into the steep-spectrum category, ~10% of the total 
non-"Ex" population. This percentage of steep-spectrum 
sources is much higher than in previous surveys, e.g. 0.5% 
in ? and 2.7% in ?. This difference is probably due to the 
increased sensitivity of both the present survey and the 
higher frequency surveys used for comparison than those 
used in these previous studies. 

Flat Spectrum Sources Flat spectrum sources are 
most likely dominated by a combination of FRI/FRII ra- 
dio galaxies both are radio galaxies whose radio emis- 
sion is powered by a central AGN, but the radio emission 
of an FRI is core-dominated and in general is less lumi- 
nous that an FRII radio galaxy, whose radio omission is 
lobe-dominated (?) - and more exotic extragalactic ob- 
jects such as BL Lacs. HII regions, which often have a 
spectral index in this frequency range of a > —0.1, SNRs, 
whose spectral indices in general fall between a ~ — 1 and 
a w —0.3 (? ?; ? ?), and PWNe, which have a radio 
spectral index of a « —0.1 (?), in M31 or the Milky Way, 
would also fall in this category. 

Sources with Spectral Curvature Sources with 
spectral curvature could be time-variable flat spectrum 
radio galaxies^, oi if tp > 0, normal galaxies (i.e galax- 
ies whose emission is not dominated by a central AGN; ? 
?), as well as AGN, where spectral turnover is most likely 
due to synchrotron self-absorption (?). In addition, if lo- 
cated in the optical disk of M31, (p > could indicate a 
foreground HII region along the line of sight, or the result 
of internal free-free absorption in a HII region or SNR. An 
extragalactic source with ip < 0, in addition to possibly 
being variable, could indicate the presence of a compact 
steep spectrum core which dominates the emission at low 
frequencies (? ?; ? ?). If the source is located in optical 
disk of M31, if < may be the result of a HII/SNR com- 
plex - if a SNR and HII region fall along the; same line of 
sight, the HII region could dominate at higher- frequencies 
and the SNR, due to its steeper radio spectrum, would 
dominate at lower-frequencies, or a diffuse particle accel- 
eration, which has been observed in young SNRs (?). 

The different kinds of sources in each of these spectral clas- 
sifications can be differentiated from each other by their 
morphological and multi-wavelength properties discussed 
in Sections 2.2 and 2.3. The criteria used to differentiate 
them is discussed in Section 3. 



2.2. Morphological Classifications 

Classifying a source based on its morphology is the 
least ambigiious method for determining the nature of the 
source used in this paper. Of the 405 sources detected, 
only the 51 "Complex" sources, shown in Figure 2, were 
sufficiently resolved to classify morphologically. All but a 
few have the radio morphology of a FRI, of a FRII, or of 
a radio jet. The classification of these sources is given in 
Section 3, and sources with anomalous morphologies are 
discussed further in Section 3.4.2. 

2.3. MulU-wavelength Classifications 

In order to classify sources using their properties at 
other wavelengths, the GLG source list was compared to 
source lists derived from targeted near infrared (IR), far 
infrared (FIR), optical, and X-ray observations of M31, 
as well as all-sky surveys in these wavebands. We looked 
for counterparts across the electromagnetic spectrum be- 
cause the different wavebands are sensitive to different 
kinds of objects, many of which could appear in the GLG 
source list. The result of these comparisons is summa- 
rized in Table ??. In this table, a counterpart is where a 
GLG source had at least one counterpart in this catalog 
and a good counterpart is the where the GLG source had 
only one counterpart in this catalog and no other GLG 
source had this catalog source as a counterpart. For "C" 
sources, which were resolved in our survey, a good coun- 
terpart could be spurious due to the large search radius 
used. For these sources, we checked if the counterpart cor- 
responded to a region of high radio emission. The number 
of false counterparts and false good counterparts were cal- 
culated by shifting the RA and DEC of every source in 
the catalog and then determining the number of counter- 
parts and "good" counterparts. In Table ??, the number 
of false counterparts and false "good" counterparts given 
is the average number of counterparts/" good" counterpart 
in the different shifts, and the error is the standard devi- 
ation. The amount of the shift was typically ±l'in both 
RA and DEC, but was larger for some catalogs (which we 
note when we discuss them). 

2.3.1. Near- Infrared 

In order to determine if a non-extended GLG source is 

an IR source, we searched for counterparts in the 2MASS 
Extended Source Catalog (ESC; ? 1)^ and the 2MASS 
Point Source Catalog (PSC; ? ?). We did not look for 
IR coimterparts of "Ex" sources because their large size, 
coupled with the high density of 2MASS sources in the 
direction of M31, resulted in a high coincidence rate, and 
because we believe that most "Ex" sources are spurious. 
We did include "ExG" sources in this analysis. For the 
ESC, the limiting distance for a counterpart between a 
PSC source and an "U" or "El" GLG source is defined to 
be: 

r = Max{3 x y(aGiG^2r^"(^Gi:G^ ^ 



3 X 



(3) 



^ If the brightness of a radio source substantially changed between the higher frequency radio surveys and our 325 MHz observation, it would 
manifest itself as curvature in the composite radio spectrum. 

* Only ESC sources within 1° of the center of M31 were included in this search. 
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where cr^^*^ is the error in RA of the GLG source, crf'^'-^ 
is the error in DEC of the GLG source, Off'-' is the major 
axis of the error clHpsc of the ESC source, and 9^^^^ is 
the minor axis of the error eUipse of the ESC source. The 
errors on position were used for the "U" and "El" GLG 
sources instead of its extent because their extents are not 
well known. For a complex GLG source, the limiting dis- 
tance for a ESC counterpart was: 

nGLG 



parameter, is -0.56, as determined using the latest cosmo- 
logical results - for the redshift of the galaxy. Using this 
redshift, we calculated the expected K-magnitude of such 
an FRI/FRII using the K-z relationship for radio galaxies 
in the 7C Redshift survey (?): 



ii: = 17.37 + 4.53 log ^-0.31(log zf 



(9) 



r = Max 



iGLG , 
'M 



(4) 

where 0^"'^ is the major axis of the GLG source, is 
the minor axis of the GLG source, and Off''' and Off'-' are 
the same as before. We used the measured size of the "C" 
GLG sources instead of the error on their positions because 
they are resolved and the IR emission need not be centered 
on the radio emission. Once this list was compiled, we re- 
jected associated we considered to be unlikely due to the 
large offset between the radio and ESC source. The re- 
maining candidates are listed in Table ??, along with the 
IR properties of the ESC counterparts. The number of 
likely associations is consistent with that expected from 
the analysis presented in Table ??. 

For the 2MASS PSC, counterparts were detected using 
the same method described above. Again, we rejected as- 
sociations where there was a large offset between the radio 
and IR source and, for resolved radio sources, no corre- 
spondence between the IR and radio emission. As seen in 
Table ??, the number of false good counterparts for "U" 
or "El" GLG sources is small. The results of this compar- 
ison, and the IR properties of the PSC counterparts, are 
listed in Table ??. 

The presence of an IR counterpart in the PSC allowed us 
to determine if a GLG source was either an FRI or FRII, 
background "normal" galaxy, a star in the Milky Way, or 
none of the above. To determine if the observed IR mag- 
nitude was consistent with that of a background FRI or 
FRII, we first calculated the redshift of the GLG source if 
it is an FRI or an FRII by simultaneously solving the fol- 
lowing equations - the luminosity-redshift (L-z) relation 
for FRI and FRIIs^ (?): 

log Lfof^ = (1.03 ± 0.22) log z + (1.92 ± 0.13) (5) 

log Lfot" = (26.97 ± 0.32) log z + (28.38 ± 0.12) (6) 

where Lf^^^ is the luminosity in Watts of an FR I galaxy 
bX y = 408 MHz, Lfof is the luminosity in Watts of an 
FRII galaxy at u = 408 MHz, and z is the redshift of the 
FR I/FR II; the luminosity-distance relationship: 

^408 = ^408 X 47rd| (7) 
where 5408 is determined by using the spectral model de- 
fined in Equation 4 in Paper I; and (II the luminosity dis- 
tance, is defined as: 

1 



As a result of this analysis, we concluded that none of 
these sources satisfied the requirement for an FRI or FRII 
radio galaxy since the required redshift were ridiculously 
high {z > 10). 

We used the IR colors of PSC counterparts to determine 
if a GLG source is a background normal galaxy, a star in 
M31, or a star in the Milky Way. Since these classes of 
objects have very different IR colors, once we determine 
the IR colors of these objects in the field of M31 we can 
classify the nature of a GLG source based on the IR col- 
ors of its PSC counterpart. The IR colors of background 
normal galaxies was determined by looking at the IR col- 
ors of all 2MASS ESC sources within 1° of M31 (top-left 
panel of Figure 3), since that catalog is believed to be com- 
posed almost entirely (^97%) of galaxies^. To distinguish 
between IR sources in the Milky Way and IR sources in 
M31, we looked at the IR color diagram of PSC sources 
within 5' of the center of M31 - expected to be dominated 
by sources in M31 (bottom-left panel of Figure 3), and 
that of PSC sources >90' from the center of M31 which 
is dominated by foreground stars and background galaxies 
(bottom-right panel of Figure 3). As Figure 3 shows, these 
two groups have a very different range of IR colors. The 
band of sources in the IR color diagram of PSC sources 
>90' is the same as that seen in IR color diagrams of stars 
in the Milky Way (?), and PSC counterparts which fall in 
this range are believe to be stars in the Milky Way. The 
results of this analysis are presented in Section 3. 



2.3.2. Far- Infrared 

We looked for Far Infrared (FIR) counterpart of GLG 
sources in the Schmidtobreick et. al., 2000 ("Knots") cata- 
log of FIR knots in M31 (?), the Xu & Hclou catalog of FIR 
point sources in M31 (?), and the IRAS Point Source Cata- 
log (IRAS-PSC; ? ?), IRAS Serendipitous Source Catalog 
(IRAS-SSC; ? ?), and IRAS Faint Source Catalog (IRAS- 
FSC; ? ?). For these catalogs, matches were determined 
using the same criteria used for radio catalogs, given in 
Equation 2. As seen in Table ??, only for the IRAS-PSC 
and IRAS-FSC catalogs are the number of counterparts 
statistically significant.^ However, even in these cases the 
offsets are rather large. Therefore, we are not able to claim 
any overlap between the observed radio and FIR sources 
in this field. However, Spitzer with its improved angular 
resolution will improve this situation. 

This equation was derived using galaxies with a redshift z<0.5. This analysis assumes that the luminosity-redshift relationship of FRIs and 
FRIIs does not change at higher redshift. 

^ Prom "Explanatory Supplement to the 2MASS All Sky Data Release" (http://www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup.html) 

^ The number of false counterparts and false good counterparts for the FIR surveys were calculated using position shifts of 10', instead of 
1' used for the other surveys, because of the large reported beam size. 



2<' 



(8) 



where c is the speed of light, Hq is Hubble's Constant (as- 
sumed to be 71 km sec~^ Mpc~^), and g'o, the deceleration 
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of M31 (?), and a XMM-Newton survey of the northern- 
half of M31 (Sergey Trudolyubov, private communication). 
The method of finding counterparts in these X-ray cata- 
logs was the same as that used to find counterparts in the 
2MASS PSC survey, described in Section 2.3.1.'' From 
the analysis presented in Table ??, a majority of RASS 
counterparts and all the counterparts from HRC survey 
of the optical disk of M31 are probably coincidental. For 
the Second ROSAT FSPC survey of M31 about one-third 
of the counterparts arc coincidental, but for the Chan- 
dra ACIS-I survey of the central region of M31 and the 
XMM-Newton survey the counterparts are probably real. 
Again, for these surveys we attempted to separate the real 
counterparts from the false ones by using the offsets and - 
for the resolved radio sources - the locations of the X-ray 
emission with respect to the radio emission. Table ?? lists 
the good X-ray counterparts of GLG sources, as well as 
their Hardness Ratios (HRl and HR2) the X-ray equiv- 
alent of IR colors - where available. HRl and HR2 are 
defined as: 

HRl = :£i±^ (10) 



HR2 



H1 + H2 + S 
H2-H1 



H2 + HI 



(11) 



2.3.3. Optical 

The GLG source list was compared to several optical 
catalogs in order to determine the nature of the source: 
the Winkler & Williams SNR catalog (WW; Ben Williams, 
private communications 2002), the Magnier et al. SNR 
catalog (Mag; ? ?), the Braun & Walterbos SNR cata- 
log (BW SNR; ? ?), the Walterbos & Braun catalog of 
Ionized Nebulae (Ion; ? ?), the Ford & Jacoby catalog of 
Planetary Nebulae (FJ PN; ? ?), and the Pellet et. al. 
catalog of HII regions (Pellet; ? ?). For all of these, we 
used the same criteria that wc used for other radio cata- 
logs (Equation 2). As seen in Table ??, the counterparts 
to GLG sources in the WW SNR catalog are likely to be 
real but the rest are most likely coincidental. Table ?? 
lists the counterparts of GLG sources in this catalog. 

In addition to these "specialty" optical catalogs, we 
compared the GLG source list with the Magnier et. al. 
BVRI catalog of optical sources in the field of M31 
(MagOpt; ? ?). While a majority of the sources in 
this catalog are stars in M31, it also contains background 
galaxies and foreground (Milky Way) stars. The criteria 
for a counterpart in this catalog is the same as that used 
for the 2MASS PSC, described in Section 2.3.1. MagOpt 
coimtcrparts of complex GLG sources were also deter- 
mined the same way they were for the 2MASS PSC, as 
described in Section 2.3.1. From the analysis present in 
Table ??, a majority of the counterparts between MagOpt 
sources and "U" or "El" GLG sources are coincidental, 
but due to the small ofi^sets between our radio sources and 
these sources, as well as the lack of optical color infor- 
mation for most of these counterparts, it is difficult to 
separate the real from the false counterparts. As a result, 
for completeness we include Table ??, which lists the GLG 
sources with a counterpart in this catalog, as well as B, V, 
R, and I magnitude of the optical coimterpart. 

We also looked for optical counterparts of GLG sources 
in unpublished Ha, [SII], and R-band images of M31 (gen- 
erously provided by Ben Williams).^ These images have a 
larger FOV than the MagOpt catalog, and the Ha & [SII] 
information helps determine if a source is a SNR. Figure 
4 shows all GLG sources with R-band counterparts, Fig- 
ure 5 shows all sources with Ha counterparts and Figure 
6 shows all sources with [SII] counterparts. In the case 
of GLG203, we believe the observed Ha emission is an 
artifact of the continuum subtraction process. 

The final classification of GLG sources with optical 
counterparts is given in Section 3. 

2.3.4. X-ray 

In order to determine if a GLG source has an X-ray 
counterpart, we compared the GLG source list to source 
lists derived from the ROS AT All-Skj Survey ("RASS") 
[both the faint (?) and bright (?) source lists], the Sec- 
ond ROSAT PSPC survey of M31 (?), a Chandra survey 
of bright X-ray sources in M31 Globular Clusters (?), a 
Chandra ACIS-I survey of the central region of M31 (?), 
a Chandra HRC survey over most of the optical disk of 
M31 (?), a deep Chandra HRC survey of the central region 

* Both the Ha and [SII] images are continuum subtracted. 

^ For the ROSAT PSPC survey, the size of the X-ray source used was the apos (Column 9) value of the source given in Table 6 of ?. For the 
Chandra ACIS survey, this value is 1" (?), for the HRC survey of the optical disk it is 10" (?), for the HRC survey of the center of M31 it is 
1" (?), for the list of Chandra sources in Globular cluster it is 5" (?), and for the XMM-Newton source list it is 3". 



with S the count rate of photons with energy be- 
tween 0.1 and 0.4 KeV, HI the count rate of photons with 
=0.4-0.9 KeV, and H2 the count rate of photons with 
E^ =0.9-2.0 KeV. 

In addition to cross referencing the GLG source list with 
the aforementioned X-ray source lists, we also looked for 
counterparts in images produced by the XMM-Newton 
survey of the optical disk of M31 (image courtesy of Al- 
bert Kong) and the Chandra-ERC survey of the southern 
disk of M31 (image courtesy of Ben Wilhams). GLG025, 
GLG040, and GLG243 were discovered to have faint uncat- 
alogued counterparts in the XMM-Newton image of M31, 
while no new counterparts were found in the Chandra- 
HRC image. 

A diverse group of sources have substantial X-ray emis- 
sion, from extragalactic objects such as FRI/ FRII radio 
galaxies (e.g. ??,??) and BL-Lacs (e.g. ?) to galactic 
objects such as SNRs, pulsars and X-ray binaries. In an 
attempt to differentiate between these sources, we looked 
at the HRl vs. HR2 diagrams of these surveys/sources. 
Unfortunately, the boundaries between these groups are 
ill-defined. However, the work of ? showed that sources 
with HR2<0 are likely to be "thermal" SNRs, remnants 
whose X-ray emission is dominated by the ejecta and 
swept-up ISM. Further classification of GLG sources with 
X-ray counterparts is given in Section 3. 

3. DISCUSSION 

Based on the information described above, we have at- 
tempted to determine the nature of the radio sources de- 
tected in this observation. We placed the GLG sources 
into five broad categories: 

• Sources in M31, 
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• Sources in the Milky Way, 

• "Background" sources (sources neither in M31 nor 
in the Milky Way), 

• Anomalous sources (source with spectral, mor- 
phological, and/or niulti- wavelength characteristics 
that defy classification), and 

• Sources for which there is insufficient information 
to classify them. 

The first four categories will be discussed in more detail in 
this Section. Based on the statistical analysis discussed in 
Paper I and on population synthesis work done by Carole 
Jackson (private communication), we expect that the fifth 
category is dominated by background radio galaxies. 

3.1. M31 sources 

The discrete radio population of a normal galaxy is dom- 
inated by SNRs, HII regions, radio pulsars, and Pulsar 
Wind Nebulae (PWNe). Each of these classes of sources 
gives different and complementary pieces of information 
regarding the properties of M31 - SNRs, pulsars, and 
PWNe give information regarding the death of past mas- 
sive stars in M31 while HII regions provide information 
on the current massive star population of M31. Detailed 
analysis of these objects can also provide information on 
the ISM and cosmic ray population of M31. 

3.1.1. Supernova Remnants, Pulsar Wind Nebulae, and 

Pulsars 

Type II SNRs and pulsars are both formed by the 
same phenomena - the massive explosion (supernova) that 
marks the death of a massive star. When the star explodes, 
it produces a large shock wave that propagates through the 
ISM, sweeping up material, that produces the SNR. Most 
of the time, a rapidly spinning neutron star (pulsar) is 
produced, which in turn can generate a highly-relativistic 
particle wind that creates a Pulsar Wind Nebula (PWN) 
(?). Table ?? lists the 11 GLG sources that we believe are 
SNR, PWN, or pulsar candidates. The columns of Table 
?? are as follows: 

Column 1: Name of source 

Column 2: Run, as defined in Section 4.1 of Paper I, 
in kpc. 

Column 3: Radio Spectral index a 

Column 4: Radio Spectral Curvature parameter 

Column 5: Classification criteria, defined later in this 
section. GLGOll also satisfies Criteria jj^2 for a 
High- Frequency Variable (defined in Section 3.3.7), 
and GLG014 and GLG036 are also CESS sources 
(defined in Section 3.4.1). 

Column 6: Size of source if at the distance of M31, in 
parsecs 

Column 7: Isotropic 325 MHz Luminosity, L325, of the 
source assuming it is at the distance of M31, in units 
of mJy kpc^. 



Column 8: Integrated radio luminosity of the source 
{v = 10 MHz-10^ MHz) assuming the source is at 
the distance of M31. This number does not take 
into account any detected spectral curvature, since 
if this source is a SNR or PWN any observed spec- 
tral turnover is most likely extrinsic to the source. 

Column 9: Emission measure, EAI , of the source in 
pc cm~^, if the spectral turnover is due to free- 
free emission. The EM is only given for "spectral 
curvature" sources with Lp > Q as, defined in Section 
2.1. Only one SNR candidate is a "spectral curva- 
ture" source. None of the PWN or pulsar candi- 
dates showed any spectral curvature. 

A GLG source was classified as either a SNR, PWNe, or 
a pulsar if it was inside the optical disk of M31 (i?M3i ^ 
27 kpc), was either a "U" or "El" source, and had one of 
the following properties: 

• A good counterpart in one of the optical SNR cat- 
alogs described in Section 2.3.3. (Criterion A in 
Table ??) 

• A good counterpart in the Pellet optical HII region 
catalog described in Section 2.3.3 and has a <-0.2. 

In the Milky Way, young SNRs are often observed 
to be near HII regions. Since the progenitors of 
Type II SNRs are short-lived, they often occur in 
regions with ongoing massive star formation, which 
often have HII regions. (Criterion B in Table ??) 

• Is a "spectral turnover" source with ^p <0. A SNR 
could have <0 as a result of diffuse particle ac- 
celeration (?), or if there is a flat spectrum source, 
most likely a HII region or PWNe, projected along 
the line of sight. In this case, the SNR would 
dominate the emission observed at lower frequen- 
cies and the flat spectrum source would dominate 
the emission at higher frequencies since the SNR 
has a steeper spectrum (in general, SNRs have a 
spectral index between a w —1.0 and a w —0.3) 
than these objects. A similar effect can be seen if a 
radio pulsar which have a median spectral index 
of a « —1.8 (?) - lies along the sight. However 
it is unlikely that a pulsar in M31 would be bright 
enough to outshine a flat-spectrum source, even at 
325 MHz. (Criterion C in Table ??) 

• A good X-ray counterpart with HR2<0, indicative 
of a thermal SNR (?). (Criterion D in Table ??) A 
pulsar or PWN could also have a good X-ray coun- 
terpart with HR2<0, in the case the radio emission 
would be dominated by the PWN or pulsar while 
the X-ray emission is dominated by the thermal 
SNR. 

• A good counterpart in the ROSAT 2nd PSPC sur- 
vey which has been classified as a "SNR" on the 
basis of previous optical work. (Criterion E in Ta- 
ble ??) 
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A source that met any of these criteria and has —1.0 < 
a < —0.3 was classified as a SNR, while a source that met 
this criteria and has a > —0.3 was classified as a PWN, 
and a source that met this criteria with a < —1.0 was 
classified as a pulsar candidate. While not every SNR, 
PWNe, or pulsar in the GLG catalog will satisfy these re- 
quirements, this set of criteria should eliminate most, if 
not all, background sources. 

The distribution of these sources over the optical disk 
of M31, as shown in Figure 9, appears to be relatively 
uniform. Having identified the SNR, PWNe, and pulsar 
candidates in the GLG source list, we compared their prop- 
erties to those of known objects of each type. 

Supernova Remnants The SNR candidates we iden- 
tified in M31 are listed in Table ?? and shown in Figure 
8. For comparison purpose, we calculated the 325 MHz 
flux densities (5325) and major axis {9m) of known SNRs 
in the Milky Way (MW), LMC, SMC and M33 if at the 
distance of M31. For MW SNRs, wo used radio spectral 
information from the Green SNR catalog (?) and distance 
information from Tables 1 and 2 of ?.^° For LMC and 
SMC SNRs, we used data presented in Table la of ?, while 
the properties of M33 SNRs was obtained from Table 4 of 
?. Figure 10 shows 6*325 vs 9m for these SNRs, as well 
as for the SNR candidates in M31 identified above. The 
GLG SNR candidates fall within the locus of points de- 
fined by SNRs in these other galaxies, implying that they 
are most likely not background sources. However the sam- 
ple of known SNRs in the MW, LMC, SMC, and M33 are 
subject to a myriad of selection effects, and there are GLG 
sources not classified as SNR candidates that have similar 
5*325 and Om- From Figure 10, the SNR candidates de- 
tected in our survey - if they indeed are SNRs - appear to 
be the brightest, and therefore probably youngest, SNRs 
in M31. These SNR candidates are brighter than known 
radio SNRs in M31, e.g. ?, which were too faint to be 
detected in our survey. Higher resolution radio imaging 
as well as deeper X ray imaging are needed to verify that 
they are SNRs and to better determine their properties. 

Of the SNR candidates identified, only GLG193 shows 
spectral turnover. If this spectral turnover is due to free- 
free absorption, the implied EM towards this source is 
2.1 X 10^ pc cm~^ — more than order of magnitude higher 
than the EM of ionized gas observed inside SNRs, e.g. 
Gas A which has an EM of - 1.7 x 10'' pc cm^^ (?). This 
possibly could be due to a forcgroimd HH region, either 
in the Milky Way or in M31, since the observed EM is 
similar to that of Galactic HH regions (?). However, the 
lack of Ha emission in this region, implies that any such 
HH region is obscured, possibly by a molecular cloud. The 
amount of absorption needed to block an HII region would 
require so much dust that M31 would not be visible if the 
obscuration was due to a Milky Way source. Therefore, 
both the HII region and the obscuring molecular cloud 
need to be in M31. However, this conclusion is extremely 
preliminary and more data is needed in order to determine 
the nature of the spectral turnover observed in this source. 



Pulsar Wind Nebulae Pulsar Wind Nebulae (PWNe) 
are magnetic "bubbles" inflated by an energetic wind 
produced by the central pulsar (?). The PWN candi- 
dates we identifled in M31 are listed in Table ?? and 
shown in Figure 8. The quintessential PWN in the Milky 
Way is the Crab Nebula, with a radio luminosity Lr of 
1.8 X 10^^ ergs s~^ (?) and is a bright emitter at X- 
ray, 7-ray, and even TeV energies. Two of the three 
PWN candidates detected in the GLG survey, GLGOll 
and GLG 198, have higher radio luminosities than the Crab 
Nebula - GLGOll is ~ 70x more luminous and GLG198 
is ~ 30x more luminous. Both GLGOll and GLG068 
have X-ray counterparts; the X-ray luminosity '^-'^ i^x) of 
GLGOll is ~ 2 X 10"^^ erg s~^ (Albert Kong, private com- 
munication), and the X-ray luminosity^^ of GLG068 is 
~ 2.1 X 10^^ erg s'^ (?). The X-ray luminosity of both of 
these sources is about an order of magnitude less than that 
of the Crab Nebula, L^'^"^ ~ 2.5 x lO'"'^ erg s"^ (?). Since 
the X-ray emission of a PWN reflects the current energy 
injection from the pulsar while the radio emission reflects 
the integrated energy input over the lifetime of the sys- 
tem, the ratio of Lx to is believed to be indicative of 
the system's age though the intrinsic scatter in Lx/Lfi 
is large (?). The Lx/Lr ratio of GLGOll and GLG068 is 
~1.6 and ^21.6 respectively, well within the Lx/Lr range 
of known PWNe (?). 

GLGOll also is a High-Frequency Variable (see Sec- 
tion 3.3.7 for more information), and has a counterpart 
in the WENSS survey whose flux is 3.5(7 lower than the 
flux measured here (numbers are in Table ??), but is not 
considered a Low-Frequency Variable using the criteria 
defined in Section 3.3.7. While PWNe are extremely dy- 
namic sources, they are not known to intrinsically vary by 
such a degree. This variation might be an example of an 
"Extreme Scattering Event" (ESE), in which extragalac- 
tic (outside the Milky Way) radio sources are observed 
to vary greatly in amplitude at frequencies v < 2.7 GHz. 
ESEs are possibly related to large radio continuum loops 
in the Milky Way (? ?; ? ?). While there are no Milky 
Way Galactic loops near GLGOll, the source does lie on 
the edge of a very large SNR/super-bubble in M31. It 
is possible that the observed variability is caused by this 
structure, although more information is needed. If this 
was the case, than the intrinsic size of GLGOll must be 
very small and therefore it could not be a PWN. 



Pulsar Candidates Of the three GLG sources identi- 
fied as pulsar candidates, GLG014 and GLG036 fall un- 
der the definition of Compact Extremely Steep-Spectrum 
(CESS) sources, which are discussed in detail in Section 
3.4.1. While these sources arc steep spectrum sources lo- 
cated inside the optical disk of M31, they are most likely 
not radio pulsars in M31 because, if so, they would be sev- 
eral orders of magnitude more luminous than any known 
galactic pulsar. In addition, GLG036 and GLG205 are 
slightly extended, which one would not expect from a pul- 
sar. Additionally, GLG205 lies far off the optical disk of 

We did not use distances from Table 3 of ? because these distances were calculated using the T,—D relationship, considered to be highly 
inaccurate. 

Measured between 0.3 and 7 keV. 

Measured between 0.1 and 2 keV assuming a hydrogen column density Nh = 9 x 10^" cm~^ and photon index F = —2. 
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M31, not expected for a radfo pulsar in M31. We conclude 
that wo did not detect any pulsars in M31, not surprising 
given the sensitivity of our survey. 

3.1.2. HII Regions 

HIT regions, the regions of ionized material that sur- 
round new O or B stars, are a major component of the 
radio emission of a normal galaxy. The radio spectrum 
emission of an HII region is flat at higher radio frequencies 
(a ~ —0.1) but turns over a lower frequencies as a result 
of free- free absorption. For roughly half of all Milky Way 
HII regions this turnover occurs at Vt > 325 MHz (?), 
so an HII region might appear as a flat-spectrum source 
or a spectral turnover source in our survey. To account 
for these two possibilities, a GLG source was classifled 
as a HII region if it was inside the optical disk of M31 
(-Rm31 < 27 kpc), not a Low-Frequency Variable or High- 
Frequency Variable (as defined in Section 3.3.7), and was 
either: 

• The cormterpart of an optical HII region and a fiat 

spectrum source, 

• A spectral turnover source with a flat spectrum, as 
deflned in Section 2.1, that has not been identifled 
as a SNR, or 

• A flat spectrum source, as deflned in Section 2.1, 
with no spectral turnover {\ip\< a^) that has not 
been identified as a SNR. 

Two sources, GLG016 and GLG088 satisfy the second cri- 
teria, while none satisfy the first or third. Unfortunately, a 
fiat-spectrum radio galaxy (FSRG) would also satisfy the 
second criteria. The presence of Far-IR emission, while 
not seen in all HII regions, would indicate that a source 
is a HII region since FSRGs are not strong FIR emitters. 
However, as discussed in Section 2.3.2, it is not currently 
possible to reliably associate FIR objects with the radio 
sources detected in this survey. As a result, we did not de- 
tect any strong HII region candidates in the GLG source 
list. Since most Milky Way HII regions when placed at 
the distance of M31 are much fainter than the detection 
limit of this survey, this is not unexpected. 

3.2. Milky Way Sources 

We believe that we have found two types of Milky Way 
radio sources in the GLG catalog, radio stars and Plan- 
etary Nebulae. Radio emission has previously been de- 
tected from a wide variety of stars: from pre-main se- 
quence K- stars (?) to 0-B stars (?). The cause of the 
radio emission is equally varied - the radio emission of O- 
B stars is expected to be dominated by free-free emission 
from stellar winds (?), while radio emission from pre-main 
sequence stars is believed to be the result of chromospheric 
emission (?). In addition, radio emission has been detected 
from binary systems (?), and possibly from stars orbiting 
a microquasar (?). 

A GLG source was categorized as a radio star if it had 
a bright optical/IR counterpart or had IR colors consis- 
tent with that of a star. Three GLG sources met this 
criteria: GLG023 - which has an / « 13 counterpart in 
the MagOpt catalog while most sources in this catalog 



have / = 18 - 20, GLG097 - which has the IR colors 
of a star, and GLG116 - which has a bright PSC coun- 
terpart {K Ki 5 as opposed to iiT ~ 13 — 15). R-band 
images of these sources is shown in Figure 4. All three 
have non-thermal radio spectrum («=— 0.89, —0.92, and 
—0.79, respectively), consistent with the processes men- 
tioned above. GLG023 has an optical counterpart in the ? 
catalog of bright stars in the field of M31, which has a B 
magnitude of 15.04±0.19. This is similar to the B magni- 
tude of its counterpart in the MagOpt catalog (14.9±0.04), 
implying the source is not variable and possibly ruling out 
a binary explanation for the radio emission. The spectral 
class of this star is unknown. GLG097 is also in a catalog 
of reference stars in the direction of M31 (?), and the spec- 
tral class of this star is also unknown. GLG116 is likely 
the radio counterpart of SAO 36591, a KO star (?). Radio 
emission from KO type stars have often been associated 
with Algol Type Binaries, which often are X ray sources 
and typically have radio emission on the order of a few 
mJys at ~ 5 GHz (?). GLG116 has no X-ray counter- 
part, and has 55400 = 22 mJy, brighter than that of most 
observed Algol-type binaries but not all, e.g Algol itself 
(?)• 

In addition to detecting radio stars, we believe that 
GLG347 may be a Planetary Nebula in the MUky Way. 
GLG347 is located in the SE corner of our FOV and has 
counterparts in the 5th Cambridge Survey (5C3.152), 3rd 
Bologna Survey (B3 0043+398), and the Miyun 232 MHz 
survey (MY 004309.7-^394955.8). On the basis of these 
detections, it has previously been identified as a radio 
galaxy (?). However, as shown in Figure 11, located to 
the SE of this source is a region of increased emission not 
classified as a source by STINV. Taken together, these 
regions suggest a shell of diameter ~45" and a fiux den- 
sity 5325=18.5 mJy. This is smaller and fainter than most 
galactic SNRs (?), and at 6 « —23°, much further off the 
galactic plane that most Milky Way SNRs, so we do not 
believe that is a SNR. However, this source has a diameter, 
morphology and flux density similar to that of planetary 
nebulae, and this is the basis of its classiflcation (?). There 
is an IR source in the center of the radio emission, shown 
in Figure 11, with the IR colors of a star which could be 
the source of the Planetary Nebula. 

3.3. Background Sources 

In this section, a "background" source is one neither in 
M31 or in the Milky Way. The population of background 
sources is dominated by galaxies whose radio emission is 
powered by a central AGN (primarily FRI and FRII ra- 
dio galaxies) - particularly sources with 5,325 > 10 mJy, 
and galaxies whose radio emission is dominated by syn- 
chrotron emission from relativistic electrons and free-free 
emission from HII regions (normal galaxies, ? ?) - ex- 
pected to dominate at ^325 < 10 mJy (Carole Jackson, 
private communication 2002). Identifying these sources is 
crucial in determining whether or not a source is in M31 
or the Milky Way, since, as shown in Section 4 of Paper 1, 
"background" sources dominate the GLG source list. 

Even though this class of sources is not the pri- 
mary motivation nor the primary focus of this paper, 
many of these are interesting objects: we detected a BL 
Lac, a galaxy merger, several High-Prequency and Low- 
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Frequency Variable sources, and have identified a Giant 
Radio Galaxy (GRG) candidate and several High-z Radio 
Galaxy (HzRG) candidates. This section will first briefly 
describe the sources which we believe are typical back- 
ground sources (FRI, FRII, and normal galaxies) and then 
disciiss the more exotic backgroimd sources in detail. 

3.3.1. FRI & FRII Radio Galaxies and Radio Jets 

The most common source with 5325 ~ 5 mJy are galax- 
ies whose radio emission is dominated by a central super- 
massive black hole (Active Galactic Nuclei; AGN). These 
sources typically fall into three broad morphological cate- 
gories: 

FRI: A galaxy whose radio emission appears core dom- 
inated (?). A GLG source was categorized as an 
FRI if it was either: 

• resolved and the radio emission is dominated 

by a single, central component (Criterion A in 
Table ??). A resolved PWN may have such a 
morphology, but we do not expect to be able 
to resolve a PWNe in M31. 

• an unresolved or elongated source with a 
counterpart in the ROSAT 2nd PSPC survey 
categorized as a "Galaxy" source which does 
not have a good IR counterpart (Griterion B 
in Table ??). As shown in Section 2.3.1, the 
detected IR emission from GLG sources is too 
bright for the source to be an FRI/FRII ex- 
planation. 

• an unresolved or elongated source outside the 
optical disk of M31 with an X-ray counter- 
part and no IR counterparts (Criterion C in 
Table ??). 

As mentioned in Section 2.3.1, no GLG source with 
a counterpart in the 2MASS catalogs had a K- 
magnitude consistent with that of FRI. Table ?? 
lists the nine GLG sources determined to be FRIs. 
One of these is also a High Frequency Variable, as 
defined in Section 3.3.7, and two are also High-z Ra- 
dio Galaxy candidates, as defined in Section 3.3.3. 
The IR counterpart of GLG122 is a star, and most 
likely unrelated to the radio source. All but two of 
these sources have X-ray counterparts. This Table 
is by no means complete since most likely a vast 
majority of unresolved sources detected outside of 
the optical disk of M31 are FRIs, it only lists those 
sources for which there is additional evidence that 
they are FRIs. 

FRII: A galaxy whose radio emission is dominated by 
radio jets emitted from the central black hole, and 
often more luminous than FRIs (?). FRIIs may 
or may not have a "core" component that corre- 
sponds to the central AGN. A GLG source was 
categorized as an FRII if it was resolved and had 
the morphology described above. Table ?? lists the 
twenty-three GLG sources determined to be FRIIs. 
GLG015 and GLG078 are also classified as anoma- 
lous sources because of Ha/IR emission detected 
near these sources. GLG078, GLG129, GLG187, 



and GLG296 are also HzRG candidates (Section 
3.3.3), and GLG187, GLG296, and GLG358 are 
Low-Frequency Variables (Section 3.3.7). GLG220 
and GLG212 form a radio triple, shown in Figure 
12. Seven FRIIs have 2MASS coimterparts, and of 
these four have the IR colors of a normal galaxy. For 
these four sources we believe that the radio emission 
is from the central AGN but the IR emission is from 
the rest of the host galaxy. The IR counterparts of 
the other three are most likely coincidental. 

Radio Jets: Radio jets are most commonly emitted 
from AGN, but tend to be larger in angular size 
than FRIIs which is why they are in a separate 
category. A GLG source was classified as hav- 
ing radio jets purely on a morphological basis. 
Figure 13 shows the five radio jets discovered in 
this survey: GLG031/GLG033, GLG045/GLG051, 
GLG054/GLG059, and GLG266/GLG269, and 
GLG270/GLG271/GLG275. Between GLG031 and 
GLG033 is an elongated 20 cm source in the Braun 
Catalog (Braun 3), not seen in our 325 MHz im- 
age, oriented perpendicular to GLG031/GLG033, 
that possibly is the host galaxy. The bright point 
sources in GLG031 and GLG033 are most likely hot 
spots in the jet (?). GLG045/GLG051 is a one- 
sided jet, where GLG045 is the central core while 
GLG051 is the jet itself. The bright point source 
in GLG051 is most likely the hot spot produced 
at the termination point of the jet (?). In the 
case of GLG054/GLG059, GLG054 has the mor- 
phology of a radio jet, and its extended emission is 
oriented towards GLG059. The GLG266/GLG269 
system also does not have the classical morphol- 
ogy of a radio jet - GLG269 is a steep spectrum 
(a ~ —2) complex source with the morphology of 
a one-sided jet, while GLG266 is a HzRG candi- 
date (Section 3.3.3) along the jet axis of GLG269. 
The GLG270/GLG271/GLG275 is a classic radio 
jet, with GLG271 most likely the central engine 
that produces GLG270 and GLG275. 

3.3.2. Normal Galaxies 

While Section 3.3.1 identified background AGN, this 
section identifies background normal galaxies - those 
whose radio emission is not dominated by a central AGN 
but from synchrotron emission emitted by relativistic elec- 
trons ( "cosmic rays" ) in SNRs and in diffuse gas interact- 
ing with the galaxy's magnetic field and free -free emission 
produced by HII regions. Normal radio galaxies have op- 
tical morphologies across the entire Hubble diagram (?). 
A GLG source was determined to be a normal galaxy if it 
satisfied at least one of the following criteria: 

Criterion I: Has an 2MASS PSC counterpart with the 
IR colors of a normal galaxy (as discussed in Section 
2.3.1). 

Criterion II: Has a counterpart in the 2MASS ESC cat- 
alog. 

Criterion III: Has a counterpart in the ROSAT 2nd 
PSPC survey that was classified as a galaxy and 
in either the 2MASS ESC or 2MASS PSC. 
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Table ?? lists the GLG sources identified as radio coun- 
terparts to normal galaxies. As shown in this table, sev- 
eral of these sources fall into other categories. In these 
cases, we believe that the observed radio and IR emission 
come from different parts of the source the radio emis- 
sion from the central AGN and the IR emission from the 
rest of the galaxy. GLG253 is located in the center of 
a z=0.293 galaxy cluster identified by the XMM-Newton 
satellite (?). 

Two of these som-ccs GLG050 and GLG203 arc 
"spectral curvature" sources as defined in Section 2.1. For 
GLG203 there is no evidence for ionized material along 
the line of sight, so the observed spectral turnover is most 
likely intrinsic. Intrinsic low-frequency spectral turnover 
has been seen in many normal galaxies, and can be ex- 
plained by a clumpy, non-thermal plasma with electron 
temperatures 200 K < Tg < 1500 K and electron density 
Tie = 0.07 — 0.10 cm"'' (?). For a typical galaxy, the opti- 
cal depth at 325 MHz of such a medium (n,, = 0.1 cm~^ 
and Te = 100 K) is ~0.002, an order of magnitude lower 
than that observed in GLG203. However, the EM of these 
galaxies is similar to that of the starburst galaxy M82 (? 
?; ? ?). The cause of the negative spectral curvature 
observed in GLG050 is unclear. 

3.3.3. High-z Radio Galaxies 

Three separate criteria, all sensitive to different kinds 
of High-z Radio Galaxies (HzRGs), were used to identify 
HzRG candidates. The first criterion used the fact that 
Ultra-Steep Spectrum (USS) sources, defined here to be 
a < —1.3, have been determined to be efficient tracers 
of HzRGs (? ?; ? ?). This is believed to be an in- 
trinsic property of HzRGs since the generic spectrum of 
a radio galaxy lobe, assumed to behave like an optically 
thin synchrotron source, steepens in spectral index by ~0.5 
around the "bend" frequency fb as one goes from low to 
high frequencies due to synchrotron losses. The observed 
i^b is proportional to (1 + z)^^, so at higher redshifts the 
bend-frequency occurs at lower frequencies, causing the 
object to have a steep spectral index when observed at 
frequencies above the redshifted bend frequency. Inverse- 
Compton scattering off the CMB also contributes to the 
steeping of the spectrum of HzRGs, since losses due to this 
process increase as (1 -|- z)^ (?). Using this information, a 
source was identified as a HzRG if it was resolved with the 
morphology of an FRI or an FRII and has a < —1.3. The 
four sources that met this "USS" criterion are listed in Ta- 
ble ??, implying a source density of ^1700 sr"^. This is 
more than 10 x higher than that measured in a comparison 
between the WENSS and NVSS surveys (?). 

The second criterion is designed to detect two sub- 
classes of potential HzRGs which have a more compact 
morphology - Gigahertz Peaked-Spectrum (GPS) sources 
and Compact Steep Spectrum (CSS) sources (?). Both 
classes are extremely powerful, compact^'"' sources with 
radio spectra that turnover at low frequencies. For GPS 
sources, this turnover occurs at Pt ~1 GHz, while for CSS 
sources this turnover generally occurs at ft ;$500 MHz (?). 
At frequencies above Vt, the spectral index of GPS sources 
is typically between a sa —0.5 and a « —0.9 (?), while 
for CSS sources it is between a w —0.5 and a ~ — 1.5 (?). 

GPS sources are <1 kpc in size while CSS generally are ~l-20 kpc 



The spectral turnover is most likely due to synchrotron 
self- absorption, but free— free absorption may also play a 
role (?). A GLG source was determined to be a CSS or a 
GPS source if it is a "spectral curvature" source with (p > 
and spectral index a <-0.5 that has not been classified as 
a "normal" galaxy (see Section 3.3.2 for details). Since 
only measurements at three frequencies are available, we 
are unable to distinguish between CSS and GPS sources, 
though a visual inspection of the spectra of these sources 
implies that most of them arc CSS sources. The sources 
that met this "CSS/GPS" criterion are listed in Table ??. 

High-resolution (Ores < 1") VLA images of CSS sources 
have shown that they are often doubles, triples, or even 
have jets associated with them (? ?; ? ?). CSS sources 
with extended emission almost always have cores with a 
fiat or inverted spectrum (—0.5 < a < 1.0) and extended 
emission with a spectral index of —1.0 < a < —0.7 (?). 
If these two components arc of comparable strength than 
it is possible that at higher frequencies (f ^1 GHz) the 
"fiat spectrum" core dominates while at lower frequen- 
cies {v ~300 MHz) the steeper-spectrum extended emis- 
sion dominates, resulting in a spectrum with (f <0. CSS 
sources of this nature were selected from the GLG sam- 
ple by requiring that the source be a "spectral curvature" 
source outside the optical disk of M31 with (p < 0. The 
sources that met this "CSS" criterion are also listed in 
Table ??. 

Table ?? lists the GLG sources that met at least one 
of these criteria. GLG296 was previously identified as a 
HzRG candidate based on its properties in the WENSS 
and NVSS catalogs (?), and several of these objects fall 
under other classifications as well. Follow up radio and 
optical observations are needed to determine if these ob- 
jects are truly High-z Radio Galaxies. 

3.3.4. Galaxy Merger 

Galaxy mergers are powerful radio sources because the 
collision triggers star formation, which in turn creates 
radio-bright sources such as HII regions and SNRs. Radio 
observations of galaxy mergers, such as the "Antennae" 
galaxies (NGC 4038/4039), have discovered that the radio 
emission from these sources is concentrated in the rcigion 
of brightest optical emission (?). GLG247 is classified as a 
galaxy merger as a result of its radio morphology, as well 
as its relationship with its IR, optical, and X-ray coun- 
terparts, shown in Figure 14. GLG247 has a relatively 
steep spectrum, a = —0.87, so it is possible that the radio 
population is dominated by SNRs like that of the "An- 
tennae" galaxies (?), but higher resolution radio images 
are needed to determine the source of the radio emission. 
GLG247 could also be an "X-shaped" radio galaxy (?). 

3.3.5. BL Lacs 

BL Lacs believed to be the radio emission from a jet 
aligned with our line-of sight emitted by a super-massive 
black hole - are some of the most powerful radio sources 
known and have been detected throughout the EM spec- 
trum, from the radio up to TeV energies (?). The X-ray 
counterpart of GLG105 was previously identified as a BL 
Lac candidate in the Einstein Slew Survey (?). In our 

size (?). 
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survey, GLG105 is a a complex source with a cometary 
inorpliology (shown in Figure 15) and a radio spectrum 
with spectral index a = —0.74 ± 0.10 and no curvature, 
but shows variability at high-frequencies (see Section 3.3.7 
for more details) . Similar variability has been observed in 
other BL Lacs (e.g. ? ?). GLG105 has optical, IR, and 
X-ray counterparts, and its broadband SED (radio to op- 
tical) is shown in Figure 15. This spectrum does not fit 
general spectral models of BL Lacs, in which a power-law 
extrapolation of the radio emission over-predicts the ob- 
served optical/IR emission (?). This discrepancy can be 
explained by variability between the different epochs of the 
observation, but is interesting nevertheless. Additionally, 
Ha emission has been detected near the "tail" of GLG105 
(shown in Figure 15), uncharacteristic of BL Lac objects. 
However, this source has a log{fx/fv) = 1-6, consistent 
with that of other known BL Lacs (?). Since GLG105 is 
located within the optical disk of M31, it is possible that it 
is a very luminous {Lx = 2.4 x 10"^* ergs s^^; ? ?) X-ray 
source in M31, but is extremely unlikely due to the rarity 
of such objects. 

3.3.6. Giant Radio Galaxies 

As mentioned in Section 3.3.1, some central AGN are 
known to produce relativistic jets which - if powerful 
enough - inflate a "cocoon" that expands first into the 
ISM of the host galaxy, and then into the surrounding In- 
tergalactic Galactic Medium (IGM). The magnetic fields 
in these cocoons produce synchrotron emission, which is 
observed in the form of radio lobes. A Giant Radio Galaxy 
(GRG) is a radio galaxy whose lobes span a projected dis- 
tance of > 1 Mpc. Since radio lobes expand with time, 
GRGs must be very old (> 10^ years), and located in 
under-dense regions relative to other, more compact, FRIs 
and FRIIs. GRGs result of their large size, very 

powerful probes of the IGM (?). 

We beheve that GLG242 and GLG260, shown in Figure 
16, located 4' away from each other, together comprise 
a Giant Radio Galaxy (GRG). (According to the search 
criteria of ?, the radio lobes need to span 5' in order for 
the source to be considered a GRG candidate. However, 
since the redshift of this object is unknown, the physical 
distance between these two objects could be > 1 Mpc.) 
GLG242 is a complex source which composed of a point 
source with an attached region of extended emission ori- 
ented in the general direction of GLG260, an unresolved 
source in our nomenclature but which, under close exami- 
nation, has extended emission in the direction of GLG242. 
GLG242 has a flux density of 5325 = 254 mJy and a 
spectral index a ~ —1.07 ± 0.42, while GLG260 is some- 
what weaker (6*325 = 176 mJy) and has a flatter spectrum 
(a = —0.81 ± 0.031). These differences arc similar to that 
seen in other GRGs (Table 3 of ? ?). When combined, 
these sources have a total flux density of 6*325 ~430 mJy 
- well within the flux range of known GRGs (Figure la 
of ? ?), and a composite spectral index 03^5° = —0.94 - 
steeper than all but a handful of known GRGs (Table A4 
of? ?). 

A more detailed analysis of these sources is hampered 
by the lack of information about their host galaxy. As seen 
in the 325 MHz and the 1.4 GHz (NVSS) radio image of 
the GLG242/GLG260 (shown in Figure 16) system, there 



is no radio source between them. While most host galaxies 
of GRGs have ^325 > 3 mJy (the 5a RMS limit in the area 
between these sources), several host galaxies of GRGs have 
325 MHz flux densities below this limit so it is not incon- 
ceivable that this is the case for GLG242/GLG260 (Table 
A3 of ? ?). In addition, there is no obvious optical/IR 
candidate for the host galaxy. There is an IR point source 
located ~20" south of the center of GLG242, but this 
source has the IR colors of a star (see Section 2.3.1 for de- 
tails), so it is probably unrelated. GLG260 has no coimter- 
parts in the Far-IR, IR, optical, or X-ray catalogs/images 
searched. A J-band image of GLG242/GLG260 (also 
shown in Figure 16) reveals the existence of a faint 2MASS 
source (2MASS:00465464-f-4033383) on the line that con- 
nects GLG242 and GLG260, but based on its IR-colors 
and its optical properties in the USNO B-1.0 catalog (?), 
it is most likely a star. 

3.3.7. Variable Sources 

Flux variability - whether at low {v < 1 GHz) or high 

{v > 1 GHz) frequencies, or on short {t^ar <1 yr) or long 
{t'uar >1 yr) timescalcs - provides a wealth of information 
about the source itself as well as on the intervening ISM. 
Studies of variable radio sources have that they typically 
vary on timescales longer than a year (?), and that source 
which vary at low frequencies {ly ^ I GHz) have differ- 
ent properties (spectral index, galactic latitude, variabil- 
ity magnitude) than those which vary at high frequencies 
(v > 1 GHz) (? ?;??;? ?). 

It is currently believed that variability at low and high 
frequencies are caused by two separate processes: 

• A process intrinsic to the source - which is almost 
always an active galactic nucleus (AGN) of some 
sort - that dominates at higher frequencies but can 
cause variability at lower frequencies. 

• A process extrinsic to the soiuce, caused by mate- 
rial in the Milky Way, that dominates variability at 
low frequencies - Interstellar Scintillation (ISS). 

This intrinsic process is believed to occur in a majority, if 
not all, flat spectrum sources - most likely QSOs or BL 
Lacs (?) - with 5i4oo > 10 mJy (?). The extrinsic pro- 
cess is believed to be Refractive Interstellar Scintillation 
(RISS), the result of radio waves from the extragalactic 
source passing through a turbulent region of the ISM (?), 
is also used to explain slow flux variations observed in 
Galactic pulsars (?). 

Using data available from the WENSS survey, we 
searched for low-frequency variables (LFVs) in the GLG 
source list. Since the WENSS source list was derived from 
observations conducted in 1991 1993 (?) while the GLG 
sources were observed in 2000, we are sensitive to LFVs 
that vary on a ~7-10 year timescale. In addition, since 
the WENSS FOV includes the entire GLG FOV, we can 
detect LFVs across the entire FOV and see if the ISM of 
M31 plays a role. A GLG source was identified as LFV if 
it satisfied one of the following two criteria: 

Criterion 1: The GLG source had a good match in the 
WENSS catalog, and had: 
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where AS'325 is the absolute value of the difference 
in 325 MHz fliix densities between the two sources, 
•^325*^ is the measured flux density of the source 
in the GLG catalog, 5';^5^^'^'^ is the measured flux 
density of the corresponding source in the WENSS 
catalog, a'§^^ is the error of S§^^ , and a^^^^^ 
is the error of S^i^^^^ 

Criterion 2: The GLG source had no counterpart in the 
WENSS catalog, and: 

Sgt^ > 30 mJy + 5 X erf (13) 

where S§^'^ and erf are defined as above. This 
criteria looks for GLG sources that should have 

been, but were not, detected in the WENSS sur- 
vey, which has a completeness limit of 30 mJy (?). 

Table ?? lists the 10 GLG sources that were classified 
as LEV, and for comparison purposes, a non-LEV GLG 
source (last line of Table). The columns are as follows: 

Column 1: GLG name 

Column 2: Morphology Type 

Column 3: Rm31, as defined in Section 4.1 in Paper I, 
in kpc. 

Column 4: Radio Spectral Index a, calculated using 

Equation 7 in Paper I. 

Column 5: Spectral curvature parameter (f, calculated 
using Equation 7 in Paper I. 

Column 6: Name of WENSS counterpart. 

Column 7: Distance between GLG source and WENSS 
counterpart in arcseconds. 

Column 8: S§.^^ , as defined above, in mJy. 

Column 9: S^f^^'^, as defined above, in mJy 

Column 10: Variation amplitude where A5325 is 

defined above and 'S^^ = Sg\^\ Sf^^^^^ /2 

Table ?? lists the 1.4 GHz properties of these LFVs. 
The columns are as follows: 

Column 1: GLG name 

Column 2: Name of 1.4 GHz counterparts 

Column 3: Distance from GLG source to 1.4 GHz coun- 
terparts, in arcseconds. 

Column 4: Flux density of 1.4 GHz counterparts in mJy. 

Column 5: 5i4oo, the average flux density of the 1.4 GHz 
counterparts, in mJy. 

Column 6: , where AS'1400 is the difference in flux 

'51400 

density of the 1.4 GHz counterparts, and 51400 is 
defined below. 

Extended sources with no radio counterparts were excluded from this 



A total of 95 GLG sourccs^^ either had a good coimter- 
part in the WENSS catalog, or had Sg\^ >30 mJy with 
no counterparts in the WENSS catalog. Of these, '--^11% 
of them are LFVs, a significantly higher percentage than 
that observed in previous LFVs searches. Additionally, 
the observed variability is higher than in previous stud- 
ies (?), especially given the high Galactic latitude of M31 
(|6| = 22°). In addition, a large number of the LFVs in 
the GLG sample are "steep" spectrum (a < —0.5) and/or 
resolved sources, different from previous LEV studies (?). 
For five of the steep-spectrum LFVs it was possible to cal- 
culate and two of them (GLG200 and GLG340) show 
no spectral curvature - unlike previously detected steep 
spectrum LFVs (?). However, it is possible that lower fre- 
quency observations of these two sources would show spec- 
tral fiattening. Of the three resolved LFVs - two of which 
are FRIIs (GLG187 and GLG296) and one is an Extended 
source (GLG009) only GLG009 and GLG296 have sufii- 
cient information at 1400 MHz to determine if they are 
high-frequency variable (HEVs) or have spectral curva- 
ture. GLG009 is a HFV, surprising because most HEVs 
are compact (?) and there doesn't appear to be a compact 
source inside GLG009. GLG296 is not a HFV, but docs 
have spectral turnover similar to that previously detected 
in LFVs. While the larger resolution of the WENSS sur- 
vey can explain a higher flux density in this survey because 
the recorded flux is inflated by extended emission in the 
region resolved out in the GLG survey, over half of the 
detected LFVs have higher fluxes in the GLG survey than 
the WENSS survey. However, the fluxes of GLG009 and 
GLG296 is lower in the GLG survey than in the WENSS 
survey so these source might not be variables. 

By using the fluxes of a GLG source's 1.4 GHz counter- 
parts, we were able to determine if a source is a HFV. The 
criteria were: 

Criterion 1: The GLG sources has two good 1400 MHz 

counterparts, and the 6*1400 of the two counterparts 
differed by more than 5crdifi , where a^is is the error 
in the difference. 

Criterion 2: The GLG sources has more than two good 
1400 MHz counterparts, and for at least one of these 

counterparts: 

1^1400 - ^1 > 5 X ^{asf + {osf (14) 

where as is the error in 6*1400 of the 1400 MHzcoun- 

terpart in question and a-^ is the error of 6*1400, the 
weighted average of all the 1400 MHz flux densities. 

Table ?? for comparison lists the six HEVs identified in 
the GLG source list, and a non-HEV GLG source (last 
line of Table). The columns are as follows: 

Column 1: Name of GLG source 

Column 2: Radio spectral index, a, as calculated in Sec- 
tion 4.3 in Paper I. 

Column 3: Spectral curvature parameter, ^p, as calcu- 
lated in Section 4.3 in Paper I. 

Column 4: Names of 1.4 GHz counterparts, 
sample since they are, most likely, false detections. 
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Column 5: Distance from 1.4 GHz counterparts to GLG 
source, in arcseconds. 

Column 6: Flux density of 1.4 GHz counterparts, 5i4oo, 
in mJy. 

Column 7 : Variation amplitude , where AS'1400 is 

■^1400 

the difference between the minimum and maximum 
51400 measurements and /S'1400 is defined above. 

Table ?? lists the 325 MHz properties of these HFVs. 
The columns are as follows: 

Column 1: Name of GLG source. 

Column 2: Morphology type. 

Column 3: Rm31, defined in Section 4.1 in Paper I, in 
kpc. 

Column 4: Name of WENSS counterpart. 

Column 5: Distance between GLG source and WENSS 
counterpart in arcseconds. 

Column 6: Flux density of source in GLG observation, 
Sg^^ in mJy. 

Column 7: Flux density of WENSS counterpart, 

S™SS in mJy. 

Colu mn 8 : Average 325 MHz flux density of source, 
5325, in mJy. 

Column 9: 325 MHz Variation amplitude, -^ais., where 

i>325 

AS'325 and 5325 are defined above. 

There are 120 sources in the GLG source list with two 
or more good 1400 MHz counterparts, and only six (5%) 
of them fit the above criteria for a HFV. This percent- 
age is slightly lower than previous searches (e.g. ? ?; ? 
?), but the amount of variation observed is similar (?). 
As mentioned earlier, only GLG009 is both a LFV and 
HFV. Two of the six HFVs were resolved - GLG009, and 
the BL Lac candidate GLG105 which was described in 
Section 3.3.5. Of the four unresolved HFVs, only one of 
them is a flat-spectrum source (GLGOll). Only one of the 
steep-spectrum HFVs (GLG004) showed any appreciable 
spectral turnover, and only at the 2cr level. 

To determine if M31 has an effect on low-frequency 
and / or high-frequency variability, we plotted ^ at both 
325 and 1400 MHz as a function of -Rmsi, seen in Figure 
17. There is no observed correlation between M31 and flux 
variation at either 325 or 1400 MHz. 

To summarize, in the GLG source list we detected 10 
LFVs and seven HFVs. A higher percentage of LFVs 
were found in the GLG survey than in other LFV searches, 
most likely the result of this survey's sensitivity to fainter 
sources. However, the variability is higher than seen in 
previous LFV surveys, and much higher than expected for 
the Galactic latitude of M31. Many of these have steep 
spectra, unlike in other surveys, and several were resolved. 
However, the HFVs detected in the GLG survey have sim- 
ilar properties seen in other HFV searches. Lastly, flux 
density variations at both 325 and 1400 MHz appear to 
be uncorrelated with structures in M31. 



3.4. "Anomalous" Sources 

In this section, we discuss sources with spectral, mor- 
phological, and/or multi- wavelength properties we do not 
understand. 

3.4.1. Compact, Extremely Steep Spectrum Sources 

Table ?? lists the 22 GLG sources which we consider to 
be Compact, Extremely Steep Spectrum (CESS) sources 
- sources that are either "U" or "El" in our radio im- 
ages, not classified as HFVs, LFVs (see Section 3.3.7 for 
details), or as HzRGs (see Section 3.3.3 for details), and 
have a spectral index a < —1.6. Images of these sources 
are shown in Figure 18. Only four of these sources have 
been detected at other wavelengths - GLG014 has an X - 
ray counterpart, GLG036 has an optical counterpart, and 
both GLG140 and GLG149 have IR counterparts, though 
these may be coincidental. This survey is not the first 
time that such objects have been detected - a comparison 
of the NVSS survey and the 365 MHz Texas survey cata- 
logs found 74 sources with a <-1.5 and four sources with 
spectral index a <-2.5 (?). This search covered 10 sr of 
sky, a much larger region than our observation, but only 
included sources with 5*365 >200 mJy (?), a much higher 
cutoff than this survey (5*325 ^3 mJy) - so the higher den- 
sity of CESS sources detected in this survey could be a 
result of our sensitivity to fainter sources. 

Two possible explanations for the observed steep- 
spectra observed are; 

Option #1: These sources have flatter spectrum than 
those calculated here, and their apparent spectrum 
is steep because of some intrinsic or extrinsic effect. 

Option #2: These sources have intrinsically steep spec- 
tra. 

These options are not mutually exclusive, and neither one 
necessarily explains all the CESS sources. 

Option #1 essentially implies that the source is vari- 
able - either intrinsically or extrinsically. As discussed 
in Section 3.3.7, such variability can either occur at low 
frequencies (in this case = 325 MHz) or at high frequen- 
cies (in this case v = lA GHz). As mentioned earlier, 
high-frequency variation, believed to be intrinsic to the 
source, occurs predominantly in "flat-spectrum" sources 
(a > —0.5; ? ?). To test this explanation, we calcu- 
lated the high-frequency variation needed to explain the 
observed steep spectrum of CESS sources assuming the 
source has an intrinsic spectral index of a = —0.5. Table 
?? shows the results of this analysis. The columns of this 
table are: 

Column 1: Name of GLG source. 

Column 2: 325 MHz flux density of GLG source in mJy. 

Column 3: Extrapolated 1.4 GHz flux density of GLG 
source, in mJy, assuming a spectral index of a = 

U.O, Oi4oo- 



Column 4: Measured (or upper limit) 1.4 GHz flux den- 
sity, in mJy, S'J^go^. For sources with counterparts 
in multiple 1400 MHz catalogs, this is the minimum 
measured value. 
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Column 5: Survey that this vahic for 5*1400 comes from. 
Column 6: ^f^-^w of this source, defined in Section 3.3.7. 

'51400 



Column 7: Average percent variation {Varuoo) of this 
source per year, equal to: 



V ar uoo = 100 x 



A5- 



1400 



til 



(15) 



omeas I 
'•^1400 I 



Rlld. tyears 



where AS'uoo is simply |S'i4oo 
the elapsed time between the GLG observation and 
1.4 GHz observation in Column 5. 



The required values for 



ASi 



5'l4 



for all but GLG046 and 



GLG152 are much higher than those observed in the previ- 
ous HFV surveys discussed in Section 3.3.7. This does not 
rule out intrinsic variability as the cause of CESS sources 
since the 1.4 GHz surveys used in the above calculation 
are much deeper than those used in previous HFV sur- 
veys. However, at 5i4oo < 1 — 10 mJy, normal galaxies are 
an increasingly dominant fraction of the radio galaxy pop- 
ulation (?), and since these sources are not expected to be 
extremely variable, it is not obvious that fainter 1.4 GHz 
sources will have higher flux variability. 

Variability at low-frequencies is also a possible explana- 
tion of the CESS population. As stated in Section 3.3.7, 
low-frequency variation is caused by two separate types of 
processes: intrinsic variation and interstellar scintillation. 
If low-frequency variation is intrinsic, then the source 
would be a Group 1 LFV (See Section 3.3.7), and pre- 
vious studies of such LFVs have found that these sources 
have "flat" spectrum (a > —0.5) (?). Therefore, we could 
analyze this possibility in a similar fashion as we did for 
the high-frequency variation possibility - using the avail- 
able 1.4 GHz data, calculate the 325 MHz of the object 
at the time of 1.4 GHz observation, and then determine 
the variability needed to explain the 325 MHz flux den- 
sity measured in the GLG survey. Table ?? presents the 
results of this analysis. The columns of this table are: 

Column 1: Name of GLG source. 

Column 2: 325 MHz flux density of GLG source in mJy. 

Column 3: Measured (or upper limit) 1.4 GHz flux den- 
sity, in mJy, S^iH' . 

Column 4: Survey that this value for 51400 comes from. 

Column 5: 325 MHz flux density of GLG source, in mJy, 
assuming a spectral index of a = —0.5, S^^. 

Column 6: of this source, defined in Section 3.3.7. 



Column 7: Average percent variation (Far325) of this 
source per year, equal to: 

A5325 



Var32h = 100 X 



til 



(16) 



f^years 
-,exp _ Qmeas\ 



where A5325 is simply IS'325 — <S'325''*| and tyears is 
the elapsed time between the GLG observation and 
1.4 GHz observation in Column 5. 



As was the case for High-Frcquency variation, the variabil- 
ity required to explain the observed steep-spectra of the 
CESS sources if they are flat spectrum, intrinsically vary- 
ing LFVs is much higher than that seen in previous LFV 
studies (e.g. ??;??) for aU CESS sources but GLG046 
and GLG152. Additionally, most Group 1 LFVs have op- 
tical counterparts, not true for most of the CESS sources. 
Sensitivity to lower 5*325 does not necessarily help this ex- 
planation, since at lower 5*325 normal galaxies dominate 
the radio population (Carole Jackson, private communica- 
tion 2002) and these sources are not expect to vary with 
such magnitude. 

If the CESS sources are LFVs, and the variation is 
not intrinsic, than the probable explanation for the low- 
frequency variability is interstellar scintillation (ISS). ISS, 
caused by variations in electron density {Sng) along the 
line of sight, typically comes in two forms: 

• Diff'ractive Interstellar Scintillation (DISS), which 
can cause flux modulations of up to 100% over 
timescales on the order of minutes and over a band- 
width of the order of a MHz (?). 

• RISS , the form of ISS believed to be behind most 
low-frequency variation (see Section 3.3.7 for more 
details), tends to cause a lower degree of flux mod- 
ulation than DISS but it occurs over a much longer 
timescale (on the order of hours) and a much wider 
bandwidth (several MHz) than DISS (?). 

Neither form of ISS is a viable explanation for the CESS 
sources. For the case of DISS, its limited bandwidth and 
short timescales make it an unlikely cause. According to 
the NE2001 model for the electron density of the Milky 
Way, the bandwidth of DISS towards M31 is ~0.4 kHz 
(?). Since the observations were done in spectral line mode 
with channel bandwidths of 10 kHz, the average effect of 
DISS over the entire bandwidth is negligible. In addition, 
the short timescales of DISS argue against it as the cause 
of the CESS sources. Though extragalactic^^ sources have 
a DISS timescale roughly 3-30 x larger than that of Galac- 
tic (Milky Way) objects (?), increasing the DISS timescale 
from minutes to <1 hour, since M31 was observed for a 
total of hour hours a 100% increase in 325 flux density over 
this timescale would be insufficient to create the observed 
steep spectra of these sources. 

RISS is unlikely to be the source of the CESS sources 
since the flux modulation typically produced by RISS is 
too small. Flux density monitoring of pulsars have; found 
that flux variation caused by RISS in distant pulsars is on 
the order of ^5%, insufficient to create the steep spectrum 
observed here. Extragalactic LFVs that vary as a result 
of RISS typically have flux variations on the order of a 
few percent (e.g. ?), again insuflacient to create the steep 
spectra observed here. 

While RISS and DISS are unlikely to be the cause of 
the CESS, this does not rule out ISS playing a role. There 
have been a few instances, often called Extreme Scattering 
Events (?), where a source increases in flux by more than 
a magnitude over a bandwidth of several MHz and over a 
timescale that ranges from minutes to hours. The most ex- 
treme example of this is PSR B0655+64, whose 325 MHz 
flux density increased by a factor of ~43 for roughly an 



In this context, "extragalactic" means outside the Milky Way, not outside both the Milky Way and M31. 
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hour over a bandwidth of ^4 MHz (?). The reason for this 
sudden amphfication is not known, but is behaved to be 
the result of ISS rather than an intrinsic process because 
this pulsar has been observed to be stable over timescale 
of minutes to weeks (?). If something similar occurred to a 
source during our observation, than the measured flux den- 
sity would be on the order of ~10x higher than the actual 
flux density. To see if this is a possible explanation for the 
CESS sources, we calculated the spectral index needed to 
explain the measured higher frequency flux densities as- 
suming the observed 325 MHz flux densities is 10 x higher 
than the intrinsic value. The results are summarized in Ta- 
ble ??. For ten of the CESS sources, the required spectral 
index is reasonable (a > —0.7). While it is a possibility, 
such extreme events as that seen in PSR B0655-I-64 are 
extremely rare - the event described above was the only 
one seen in a long-term pulsar monitoring project (?) - 
so it is unlikely that something similar happened to many 
objects during our observation. However, there is prece- 
dence for this sort of behavior in the direction of M31 - in 
the 36W survey, it was observed that the flux density of 
5C3.132 decreased by a factor of five in eight years (?). 

As a result of the analysis above, we do not believe that 
variability - whether intrinsic or the result of ISS - does 
a good job of explaining the observed steep spectrum for 
most of the CESS sources. Therefore, we are left with Op- 
tion #2, that the CESS sources (or at least some of them) 
are intrinsically steep-spectrum soiirces. The most com- 
mon steep spectrum sources are radio pulsars, High-z Ra- 
dio Galaxies (HzRGs), and Relic Radio galaxies (RRGs). 

If a CESS source is a pulsar, then it must be imresolved 
and is in either M31 or the Milky Way. To test this pos- 
sibility, we assume that a CESS source in the optical disk 
of M31 {Rm3i < 27 kpc) is a pulsar in M31 - meaning it 
has a distance of 780 kpc, while for a CESS source outside 
the optical disk of M31 we assumed it is a pulsar in the 
Milky Way - implying a distance of ~10 kpc. To deter- 
mine if the pulsar possibility is reasonable, we compared 
the luminosity of these sources with that of the brightest 
galactic pulsars - at = 400 MHz it is PSR B1302-64 with 
a beamed luminosity^^ L400 of 26100 mJy kpc^, while at 
1^ = 1.4 GHz the brightest pulsar is PSR B0736-40 with a 
luminosity ii4oo of 9700 mJy kpc^ (numbers taken from 
the ATNF Pulsar Database). The results of this analysis 
are shown in Table ??, whose columns are as follows: 

Column 1: GLG Name 

Column 2: i?M3i in kpc, as defined in Section 4.1 in Pa- 
per I. 

Column 3: Z/400 in niJy kpc^ 
Column 4: L1400 in mJy kpc^ 

Column 5: af^^, the spectral index required for L400 of 
the source to be equal to that of PSR B1302-64. 

Column 6: ct]^^^, the spectral index required for L1400 
of the source to be equal to that of PSR B0736-40. 

Both 1^400 and L1400 are the beamed luminosity. The lu- 
minosities of CESS sources which might be pulsars in M31 



are higher, by at least an order of magnitude, than that 
of the brightest observed Galactic pulsars. The spectral 
index required to make these luminosities consistent with 
that of the brightest 400 MHz pulsar are much steeper 
than that of any known pulsar, those for some of these 
sources the spectral index required to make the luminos- 
ity consistent with that of the brightest 1.4 GHz pulsar 
is similar to that of the steepest spectrum radio pulsars 
(?). For CESS sources outside the optical disk, their radio 
luminosities at a distance of 10 kpc is consistent with that 
of Milky Way pulsars in the Milky Way, but are unlikely 
to be pulsars because of their high Galactic latitude. It is 
extremely premature to suggest that any of these sources 
may be pulsars, since deep optical imaging is needed to 
rule out the possibility that a source is a HzRG/rclic ra- 
dio galaxy, objects expected to be faint optical sources (?). 

Another possibility is that the CESS sources are pulsars 
than emit giant pulses, which have been detected from a 
wide variety of pulsars from young pulsars like the Crab 
Pulsar (?) to the millisecond pulsar PSR B1937-h21 (?). 
Even though giant pulses can be as much as 4000 x brighter 
than the average pulse, they do not significantly increase 
the integrated flux density (?). The most extreme case of 
this is PSR B0540-69 in the LMC, which emits giant pulses 
but remains undetected in continuum 1400 MHz observa- 
tions. This pulsar has a spectral index aHo^ < —4.4 (?), 
a similar limit as for many of the CESS sources in Table 
??. 

The second possibility for the intrinsic steep-spectrum 
explanation is that they are High-z Radio Galaxies 
(HzRGs), which were discussed in Section 3.3.3. In that 
section, the steep-spectrum HzRG candidates discussed 
were all resolved radio sources with the morphology of an 
FR II. However, HzRG candidates often are much smaller 
than those sources - the ? survey of Ultra Steep Spec- 
trum (a < —1.3; USS) sources found that USS sources 
have a constant median angular size of ~6" between 
Si4oo=10 mJy-1 Jy (?). Therefore, we expect to detect un- 
resolved USS sources in this survey. The density of CESS 
sources in this survey (~9500 sr~^) is much higher than in 
other surveys, e.g. the density of USSs in the WENSS and 
NVSS survey is ~151 sr~^ (?). A higher density of steep- 
spectrum sources in the GLG survey is not completely un- 
expected since our survey is substantially deeper than the 
WENSS survey (a flux limit of S'325 ~3 mJy as opposed 
to 5*325 ~18 nijy) and the 1.4 GHz surveys used for com- 
parison are much deeper than the NVSS survey. However, 
as mentioned earlier, starburst galaxies dominate the pop- 
ulation of sources with 6*1400 ^10 mJy - which are not 
expected to have spectra as steep as these sources. 

The third possibility for an intrinsically steep spectrum 
is relic radio galaxies. These are objects which used to 
be FRII radio galaxies, but the central AGN turned off 
causing the rclativistic electrons insidcs the lobes of these 
galaxies to lose energy through synchrotron emission and 
inverse Compton Scattering off the Cosmic Microwave 
Background (CMB). The higher energy electrons lose en- 
ergy faster than the lower energy ones, steepening the elec- 
tron energy distribution and, consequently, the observed 
radio spectral index (?). These sources can have spec- 
tral indices as steep as— 6<a<— 5, but tend to be 



^® Since radio emission from pulsars are not emitted isotropically, we used L = S x not L = S x imdP to account for beaming. 
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rather large with Om > 10" (?). Only GLG006, GLG022, 
GLG036, GLG059, and GLG305 have 9m > 10", with 
GLG036 the largest with 0m = 16.6". There is an R- 
band source on the edge of GLG036, but it is unlikely the 
two arc related. Also, GLG036 is one of the weakest CESS 
sources detected. 

For the reasons stated above, we do not believe that is 
possible at this time to definitively state the nature of the 
observed CESS sources, and feel they are most likely a mix 
of the possibilities described above. Deeper optical imag- 
ing and another 325 MHz radio observation of M31 would 
be very helpful in solving this mystery: optical imaging 
would determine whether or not these sources are galax- 
ies (HzRGs tend to have i? < 20 mag), and another radio 
image would put a strong constraint on the variability of 
these sources. These sources are much fainter than previ- 
ous steep-spectrum samples, and therefore are potentially 
extremely interesting. 

3.4.2. Anomalous Morphology and Multi-Wavelength 
Characteristics 

In this section, we discuss sources with anomalous radio 
morphologies and/or multi- wavelength properties. 

GLG015 GLG015 is a complex source located on the 
edge of the optical disk of M31 with a spectral index of 
a = —0.99 ± 0.08 and no spectral curvature, and has the 
morphology of a compact FRII. However, GLG015 is bor- 
dered on two sides by two clumps of Ha and [SII] emis- 
sion (see Figures 5 and 6) and is orientated such that it 
fits between them without overlap. This emission has no 
counterparts in the R-band image of this region, or in 
any of the Far-IR, IR, and X-ray catalogs, but the bright 
spots in Ha/ [SII] correspond to stars in the MagOpt cat- 
alog. GLG015 most likely is a FRII coincidcntally located 
in this position and orientation, if it is connected to the 
Ha/ [SII] emission around it, the nature of this source is 
unknown. 

GLG065 GLG065 is a complex located just outside the 
optical disk of M31. GLG065 has the most puzzling mor- 
phology of all sources in the GLG source list - as shown 
in Figure 2 it appears to be a collection of four, maybe 
five points sources arranged somewhat symmetrically. It 
is unlikely that this source is a gravitational lens due to the 
the wide separation of the sources (at ~30" this would be 
one the most widely separated gravitational lenses known), 
and the lack of optical/IR emission expected from the lens- 
ing object. It does have a passing resemblance to 3C315, a 
classic "X-shaped" radio source (?). Another possible ex- 
planation is that GLG065 is a collection of SNRs in M31, 
although the lack of Ha emission in this region makes this 
explanation unlikely. 

GLG078 GLG078 is a complex, steep spectrum (a = 
—1.57) source with significant spectral turnover (ip = 0.10) 
just beyond the optical disk of M31 and the morphology 
of an asymmetric FRII. There is a 2MASS PSC source lo- 
cated next to the center of the radio emission - but there 
is no radio emission from the specific location of the IR 
source, and it has the IR colors of a star in the Milky 



Way. The R band image of GLG078 (shown in Figure 4) 
shows the IR source mentioned above as well as a source 
located at the SE edge of the radio emission. We believe 
that GLG078 is an FR II and that the R-band source is 
coincidental. Higher-resolution radio images are required 
to better understand the nature of this source. 

4. CONCLUSION 

In this paper we have presented the results of a four 
hour 325 MHz survey of M31. A statistical analysis of the 

405 radio sources detected in this survey, done in Paper 
I (?), revealed that most of these sources are background 
radio galaxies. In this paper we present our attempt to de- 
termine the nature of these source using their observed ra- 
dio spectral properties, radio morphologies, and their Far- 
IR/IR/Optical/X ray properties of the detected sources. 
With this information we classified 112 out of the 405 de- 
tected sources, as summarized in Table ??. 

Sources in the GLG catalog could be placed into four 
broad categories: sources in M31, sources in the Milky 
Way, extragalactic sources, and anomalous sources. In 
M31, we have identified five supernova remnant candidates 
and three pulsar wind nebula candidates, given in Section 
3.1. The supernova remnant candidates have properties 
(extent and radio luminosity) similar to that of known su- 
pernova remnants in the Milky Way, M33, LMC, and SMC 
as shown in Figure 10. Two of the pulsar wind nebula can- 
didates are more luminous than the Crab Nebula, but only 
by a factor of a few. Further observations are needed to 
confirm the classification of these sources and determine 
the morphology and basic physical parameters. The ob- 
jects are the brightest of their class in M31, and as such 
are individually interesting. 

In the Milky Way, we believed we have identified three 
radio stars based on optical and near-IR information and a 
Planetary Nebula based on radio morphology, as described 
in Section 3.2. A wide variety of extragalactic sources were 
identified, ranging from the ordinary (FRI and FRII radio 
galaxies, star-forming galaxies, and radio jets) to the ex- 
otic (a BL Lac candidate, a galaxy merger, a Giant Radio 
Galaxy candidate, and High-z Radio Galaxy Candidates, 
including Compact Steep Spectrum and Gigahertz Peaked 
Spectrum radio galaxies), using the criteria described in 
Section 3.3. Last but not least, we have identified many 
anomalous sources - ranging from highly variable sources 
(Section 3.3.7) to compact {6m ^ 10") extremely steep 
spectrum (a < —1.6) sources (Section 3.4.1) to source with 
anomalous morphologies and/or multi- wavelength proper- 
ties (Section 3.4.2). 

To summarize, in just a four-hour 325 MHz observation 
of M31 we have potentially identified the brightest super- 
nova remnants and pulsar wind nebulae in these galaxy, as 
well as identified a class of steep spectrum sources whose 
nature is unknown. In this paper, we have a detailed a 
technique for separating background radio sources from 
those intrinsic to M31, which will be extremely useful in 
future radio surveys of this galaxy. While much remains 
unknown about the radio population of M31, this survey is 
an important step in understanding the properties of this 
galaxy. 
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Fig. 1. — The location of stccp-spcctrum sources (upper left), flat- 
spectrum sources (upper right), and sources with spectral turnover 
(bottom), as defined in Section 2.1. The names are to the right of 
the source location, and the optical image of M31 is from the Palo- 
mar All Sky Survey. As mentioned in Section 2.1, the clustering of 
sources with spectral turnover around the optical disk of M31 is a 
selection effect, further described in the text. 

Fig. 2. — Greyscale images of all Complex GLG sources overlaid 
with contours. Ellipses denote the size and orientation of the source 
as determined by the MIRIAD task SFIND. All the GLG sources are 
roughly the same size because SFIND determined the properties us- 
ing the smoothed maps described in Section 2.2.2 in Paper I. The 
intensity scale is in Jy beam^^. The values of the contours change 
for each source, and correspond to 4 1 where A is the peak flux of the 
source, and i is a counter that begins at i = 1 and increases by 1 for 
sources with — i — < 10, by 2 for sources with 10 < — - — < 20, 

"RMS ~ CTRMS ~ 

and by 4 for sources with — - — > 20, until 4^ < auMS- 
■' <^RMS » mvio 

Fig. 3.— The IR colors of all ESC sources within 1° of M31 
{top-left), of good PSC counterparts of GLG sources (top-right), 
of PSC sources within 5' of M31 (bottom-left), and of PSC sources 
>90' away from M31 (bottom-right). In the top— left graph, dia- 
monds are good ESC counterparts of GLG sources and crosses are 
good PSC counterparts of GLG sources, and in the bottom graphs 
crosses indicate good counterparts of GLG sources. 

Fig. 4. — R-band (Greyscale) image of GLG sources with op- 
tical counterparts, overlaid with radio contours (only for Complex 
sources) and ellipses that show the extent and orientation of the 
GLG source. The colorbar is in ADU counts. Since the R-band im- 
ages were not calibrated, the conversion between ADU counts and 
flux is unknown. (Ben Williams, private communication) The con- 
tours change for each source, and correspond to ^ , where A is the 
peak flux of the source, and i is a counter that begins at i = 1 and 
increases by 0.5 for sources with — - — < 10, by 2 for sources with 

— ^ — < 20, by 3 for sources with 20 < — ^ — < 50, and by 10 for 

"RMS ~ "RMS ~ 

sources with 50 < — ^ — , until 4- < aaMS- 

"RMS' ' rtiuo 

Fig. 5.-- Ha image of GLG sources with an Ha counter- 
part. The colorbar is in ADU counts, and the conversion to flux 
is 3.28 X 10-14 "^j^y ^ (Ben Williams, private communication). 
The contours change for each source, but follow the formula given 
in Figure 4. 

Fig. 6. — [SII] image of GLG sources with a [SII] counter- 
part. The colorbar is in ADU counts, and the conversion to flux 
is 2.89 X 10~-^^ "^^^ ^ (Ben Williams, private communication). 
The contours change for each source, but follow the formula given 
in Figure 4. 

Fig. 7. — Radio contour images of all the Complex sources with 
good X-ray counterparts. The cross represents a Chandra HRC 
sources, while the ellipses represent ROSAT PSPC sources. The 
contour values change for each source, and correspond to y, where 
A is the peak, flux of the source, and i is a counter that begins at t = 1 
and increases by 1 for sources with — i — < 10 until 4- < "Rms 

■' "RMS — I 2 ' 

and by 2 for sources with 10 < — ^ — < 30, and by 4 for sources 

"RMS — 

with — - — > 30, until 4 < cuMS- 
"RMS ' I nivjo 



Fig. 8. — 325 MHz Greyscale images of SNR (top two rows) and 
PWN (bottom row) candidates in the GLG source list. The inten- 
sity scale is in Janskys. The ellipse corresponds to the orientation 
and extent of the GLG source. 

Fig. 9.— The position of M31 SNR and PWN candidates with 
respect to the optical disk of M31. Optical image courtesy of the 
Palomar All-Sky Survey. 

Fig. 10.— 5326 vs 6»M of Milky Way (Xs), M33 (diamonds), LMC 
(triangles), and SMC (squares) supernova remnants if placed at the 
distance of M31. SNR candidates in the GLG sources are denoted 
by the thick "-|-" signs. Arrows indicate the upper limit on 6m for 
GLG SNR candidates. The line represents the observational limits 
of the GLG survey - only SNRs with a flux density 5325 > 3 mjy 
beam~^ will be detected. 

Fig. 11. — 325 MHz Radio image (left) and J-Band image (right) 
of GLG347. The black contour indicates the location and extent of 
GLG347 and the grey contours correspond to ^325 = 1,2,3,4 and 5 
m Jy beam" ^ . 

Fig. 12.— 325 MHz grey scale of Radio Triple GLG220 and 
GLG212. The black ellipses indicate the location and extent of 
the GLG sources, and the white contours correspond to 5325 = 
2, 3, 4, 5, 7 mJy heam~^. 

Fig. 13.— Radio Jets GLG031/GLG033 (upper-left, 
previous page), GLG045/GLG051(Mppe7 — right, previous page), 
GLG054/GLG059 (bottom-left, previous page), GLG266/GLG269 
(bottom-right, previous page), GLG270/GLG271/GLG275 (this 
page). The black ellipses are the GLG sources plus IBraun 3 (middle 
ellipse in GLG031/GLG033), and the white contours are the same 
as in Figure 12. 

Fig. 14. — Radio (left) and Palomar All Sky Survey optical image 
(right) of GLG247 overlaid with radio contours (white) whose levels 
are S32S = 2,3,4,5,7 mJy beam~^. The black ellipse corresponds 
to the size and orientation of the GLG source. 

Fig. 15. — Radio image of BL Lac candidate GLG105 (uppei — left), 
Ha image of GLG105 (uppei — right; courtesy of Ben Williams), and 
broadband SED of GLG105 (bottom). The black contours in the 
top two images are 5*325 = 2.5, 5, 15 mJy beam^-'^. The arrows in 
the bottom image arc upper-limits, while the diamond and triangle 
are the expected value for an FRI and FRII, respectively, using the 
procedure described in Section 2.3.1. 

Fig. 16. - NVSS image (top-right) , Palomar All-Sky Survey image 
(Lop-leJL). and 325 MHz radio image of Giant Radio Galaxy candidate 
GLG242/GLG260 (bottom). White contours are 5325 = 2.5,5, ...,15 
mJy beam^ ^ . 

Fig. 17. — Raw (top-row) and Binned (bottom-row) Distribution 
of ^ vs. Rm31 at 325 MHz (left) and 1400 MHz (right). 

Fig. 18. — 325 MHz images of CESS sources. The ellipses corre- 
spond to the size and orientation of the GLG source. 
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